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8 1 Abstract 

Q 

We review the theoretical status of the lifetime ratios T B +/rB d and tbJtb a and of 
the mixing quantities AM S , AT S and <p s . We show that the ratio AY S /AM S can be 
determined with almost no non-perturbative uncertainties. Finally we explain how 
this precise determination of the standard model values can be used to find possible 
new physics contributions in AM S , AT S , AT S /AM S and a s j s . Combining the latest 
experimental bounds on these quantities one already gets some hints for new physics 
contributions. 

o 

^ '. 2 Introduction 

a: 

o ' Inclusive decays (see e.g. [1] or [3] and references therein) and lifetimes of heavy 

mesons can be calculated within the framework of the so-called heavy quark expansion 
(HQE) [3111]. In this approach the decay rate is calculated in an expansion in inverse 
powers of the heavy b-quark mass. 

H A 2 A 3 

P.. r = r + — r 2 + — r 3 + . . . (1) 

mi ml 

T represents the decay of a free heavy b-quark, according to this contribution all 
b-mesons have the same lifetime. The first correction arises at order 1/mf, they are 
due to the kinetic and the chromomagnetic operator. At order \jm\ the spectator 
quark gets involved in the weak annihilation and Pauli interference diagrams [3j [5]. 
This contributions are numerically enhanced by a phase space factor of 167r 2 . Each 
of the Tj contains perturbatively calculable Wilson coefficients and non-perturbative 
parameters, like decay constants or bag parameters. Unfortunately the theoretical 
predictions for the decay constants vary over a wide range: quenched lattice de- 
terminations for fs a tend to give values of (9(200) MeV, while recent unquenched 
calculations with 2+1 dynamical light flavors give values around 260 MeV - for a 
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more detailed discussion see [HI [71 [8] . Since lifetime differences depend quadratically 
on the decay constants, going from 200 MeV to 260 MeV results in an increase of 
70%. Here clearly theoretical progress is necessary to pin down the error on the decay 
constants considerably 

In view of several new theoretical and experimental developments we update the 
numbers present in the literature (see e.g. [9]). 



3 Lifetimes 

The lifetime ratio of two heavy mesons can be written as 

!l = i + 4(rf + fir<" + ...) + ^(rr + ^r<" + . ..)+... (2) 

t 2 m b \ 4:71 J ml \ 47r / 

If one neglects small isospin or SU(3) violating effects one has no 1/ml corrections 
and a deviation of the lifetime ratio from one starts at order l/m|. 



3.1 TB+/TB d 

The leading term I3 has been determined in |3j[T0]. For a quantitative treatment of 
the lifetime ratios NLO QCD corrections are mandatory - 13 has been determined in 
[TTl [T2] . Subleading effects of 0{l/m b ) turned out to be negligible [13]. The matrix 
elements of the arising four-quark operators have been determined in [T4]. Using the 
result from [TT] 



t(B- 



0.0325 



t(B+) 



cb\ 



T(B + ] 
m b 



(3) 



B 



1.653 ps 10.04/ V4.8GeV7 \ 200MeV 



1.0 ±0.2) B x + (0.1±0.1)5 2 - (18.4 ± 0.9) ei + (4.0 ± 0.2) e 2 J + 6 1/m 

one gets with the matrix elements from Becirevic [2] (Bi = 1.10 ± 0.20; B 2 = 
0.79±0.10; ei = -0.02±0.02; e 2 = 0.03±0.01) and the values V cb = 0.0415, m b = 4.63 
GeV and f B = 216 MeV g5\: 



r(B+ 



1.063 ±0.027. 



J NLO 



which is in excellent agreement with the experimental number [T6J [T7j 

~t(B+ 



r(B° d ) 



1.071 ±0.009. 



(4) 



(5) 



1 In the case of T\ b Jt s d these effects are expected to be of the order of 5%. 
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From Eq. (j3J) one clearly sees that a precise knowledge of the color octet bag pa- 
rameters £1 and €2 - these parameters are of order 1/N C - is mandatory since their 
coefficients are numerically enhanced. Here clearly more work has to be done. 

3.2 T Bs /T Bd 

In the lifetime ratio r^ s /r^ d a cancellation of weak annihilation contributions arises, 
that differ only by small SU (3)- violation effects. One expects a number that is close 

to one PH H21 esi nsi 

4^4 = 1.00 ±0.01. (6) 

T { B d) 

This expectation is confirmed by experiment [201 E] 

4|H = 0.957 ±0.027, (7) 

T \ B d) 

although more precise experimental numbers would be very desirable. 

3.3 t b + 

The lifetime of the doubly heavy meson B c has been investigated in [21] 

r(B c ) = 0.52±g;i|ps. (8) 

In addition to the b-quark now also the c-charm quark can decay, giving rise to the 
biggest contribution to the total decay rate. The current experimental number [22] 

t{B c ) = 0.469 ±0.027ps (9) 

agrees nicely with the theoretical prediction, but it has much smaller errors. Here 
clearly some theoretical improvements are necessary to pin down the error. 



4 Mixing Parameters 

In this section we briefly investigate the status of the mixing parameters. For a more 
detailed review we refer the interested reader to [TJ. 

The mixing of the neutral B-mesons is described by the off diagonal elements Ti2 and 
M12 of the mixing matrix. Ti2 stems from the absorptive part of the box diagrams - 
only internal up and charm quarks contribute, while M\2 stems from the dispersive 
part of the box diagram, therefore being sensitive to heavy internal particles like 
the top quark or heavy new physics particles. By diagonalizing the mixing matrix we 
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obtain the physical eigenstates B B and Bl with defined masses {Mr, Ml) and defined 
decay rates (T#, T^) in terms of the flavor eigenstates B s = (bs) and B s = (bs): 

B H :=pB + qB, B L :=pB-qB with \p\ 2 + \q\ 2 = 1 . (10) 

The calculable quantities |M 12 |, |r 12 | and = arg(— M 12 /r 12 ) can be related to three 
observables: 

• Mass difference: 



AM := M H -M L = 2|M 12 | (\ + \j^ 2 ^ <t> + -J (H) 

|M 12 | is due to heavy internal particles in the boxdiagrams like the top-quark 
or SUSY-particles. 

Decay rate difference: 

Ar := T L -T H = 2|r 12 | cos0 ^1 - sin 2 <f>+..^ (12) 

1 17*12 1 is due to light internal particles: particles, like the up- and the charm- 
quark. It is therefore very insensitive to new physics contributions. 

Flavor specific/semileptonic CP asymmetries: 

A decay B q — > f is called flavor specific, if the decays B q — > / and B q ^> f are 
forbidden and if no direct CP violation occurs, i.e. \(f\B q ) \ = \ (f\B q )\. Some 
examples are B s — ► D~7T + or B q — > Xlv (therefore the name semileptonic CP 
asymmetry). The flavor specific CP asymmetry is defined as 

g 

P 



r s ff (f) - /) - r s 9 (f) - / 

a fs ~- 



T[B q (t)^f)+T[B q (t)^f 

= 1m— — = — — tan d> . 
My, AM 

4.1 Mass difference 

Calculating the box diagram with internal top quarks one obtains 

M 12 , q = ^(V*V tb ) 2 M^S (x t )B B J 2 B M Bq f, B (13) 

The Inami-Lim function So(x t = fnf/M^) [23J is the result of the box diagram without 
any gluon corrections. The NLO QCD correction is parameterized by r) B w 0.84 
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[21]. The non-perturbative matrix element of the operator Q = (bq)v-A(bq)v~A is 
parameterized by the bag parameter B and the decay constant /# 

(B q \Q\B q } = ^B Bq f B M Bq . (14) 

Using the conservative estimate fs s = 240 ± 40 MeV [7] and the bag parameter B 
from JLQCD (25] we obtain 

AM, = 19.3 ±6.4 ± 1.9 ps" 1 (15) 

The first error stems from the uncertainty in fs s and the second error summarizes 
the remaining theoretical uncertainties. The determination of AM^ is affected by 
even larger uncertainties because here one has to extrapolate to the small mass of the 
down-quark. The ratio AM s /AMd is theoretically better under control since in the 
ratio of the non-perturbative parameters many systematic errors cancel. 

This year also AM S was measured, leading to the pleasant situation of having 
very precise experimental numbers at hand J2D1 [251 12~T] 

AM d = 0.507 ± 0.004 ps^ 1 , (16) 
A M s = 17.77 ± 0.12 ps" 1 . (17) 

To be able to distinguish possible new physics contributions to AM S from QCD 
uncertainties much more precise numbers for fg a are needed. 



4.2 Decay rate difference and flavor specific CP asymmetries 

In order to determine the decay rate difference of the neutral B-mesons and flavor 
specific CP asymmetries a precise determination of T\2 is needed. With the help of 
the HQE T 12 can be written as 

r» = 4 (if + ?rf> + ...) + * (rf > + firf + ...) + ... (is) 

ml V 47r / ml \ 4tt J 

The leading term I3 was determined in [28]. The numerical and conceptual im- 
portant NLO-QCD corrections (I3 ) were determined in [221 [301 ED]- Subleading 
1/m-corrections, i.e. I4 were calculated in [191 [32] an( i even the Wilson coefficients 
of the l/m 2 -corrections (r^) were calculated and found to be small |33j . 

Besides the already known operator Q in the calculation of T12 three additional 
operators arise Q, Qs and Qs- The tilde stands for a color rearrangement and index 
S corresponds to a S-P Dirac structure instead of the V-A-structure. Q = Q and it 
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can be shown [T9J, [221 CTJ that a certain combination of Q, Qs and Qs is suppressed 
by powers of l/m b and a s 

Q = Q and R = Qs + (*iQs + ?<2 = O ( — , a.) (19) 

with ai = 1 + O (a s ), for more details see [7j. In the literature [28l fl9l [29| |30l |3T] 
always was eliminated - with the help of Eq. ( [191) - and one was left with the 
operator basis {Q,Qs}, which we call in the following the old basis. Working in the 
old basis one finds several serious drawbacks: 

• An almost complete cancellation of the coefficient of the operator Q takes place, 
while the operator is dominant. So in the ratio Ar s / AM S the only coefficient 
that is free of non-perturbative uncertainties is numerically negligible. 

• The 1/m corrections are abnormally large - all contributions have the same sign. 

• The a s -corrections and the remaining /i-dependence is unexpectedly large. 

In [7J it was found, that expressing T 12 in terms of the new basis {Q, Qs} one gets a 
result, that is free of the above shortcomings. The change of the basis corresponds 
to throwing away certain contributions of 0(a 2 ) and 0(a s /m b ), which is beyond the 
calculated accuracy. For our new determination of Ti2 we also use the MS-scheme 
|34j . besides the pole scheme for the b-quark mass. Moreover we sum up logarithms 
of the form zlnz - with z = m 2 c /m1 - to all orders, following [11] and of course we 
have to include also sub leading CKM-structures to determine a/ s , as done in 
In the old basis one obtains 



Ar, 

AM S 
with 



fB a \ 2 



240 MeV 



0.0025 + 0.0945^ - (0.0335^ + 0.0195^ + 0.0055/? 



10" 



0.9 + 40.9^f - ( 14.4-J^ + 8.5%^ + 2.1^ 
B \ B B B 



(20) 
,(21) 



K /l/f2 

(B s \Q s \B s ) = --B' s f 2 B M Ba , B' X :=B X B L (22) 
3 (m b + m s ) 



In Eq. (|2~T1) we have explicitly shown the dependence on the dominant 1/m operators 
i?2 and Rq (see [T91 [7] for the definition). The remaining power corrections are 
summarized in the coefficient of Br. One clearly sees that the cancellation in the 
coefficient of B leads to the undesirable situation, that the only coefficient in AT / AM 
that is free of non-perturbative uncertainties is negligible. 
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This changes however dramatically if one uses the new basis 



Ar 5 . 

AM, 



with 



240 MeV 



0.1055 + 0.0245^ - (0.0305^ - 0.00QB Ro + 0.0035^] 



10" 



46.2 



.B'c 



s 



+ 10.6-f 



B 
1 



'l3.2^-2.5^ + 1.2^ 
B B B 




{B S \Q S \B S } = -B' s flM Bs . 



(25) 



Now the dominant part of AT/ AM can be determined without any hadronic uncer- 
tainties! 

Using the non-perturbative parameters from [25| [35] . we obtain the following final 
numbers (see [7] for the complete list of the numerical input parameters) 



Ar, 



Ar s 

AM, 



(0.096 ±0.039) ps" 1 
(2.06 ±0.57) • 10~ 5 , 
(49.7 ±9.4) • 10~ 4 , 
0.0041 ± 0.0008 = 



Ar 

— - = AT s -r Bd =0.147 ±0.060, (26) 

J- s 

(27) 
(28) 

0.24° ±0.04. (29) 



The composition of the theoretical error of Ar is compared for the use of the old and 
the new basis in Fig. ([1]). The by far dominant error comes from the decay constant 
f Bg , followed by the uncertainty due to the power suppressed operator R 2 and the 
remaining //-dependence. In this case the theoretical improvement due to the change 
of basis is somehow limited by the huge uncertainty due to f Bs , which is the same in 
both bases. 

This changes if one looks at the composition of the theoretical error of AT / AM 
in Fig. ([2]). Since now f Bg cancels the dominant error comes from the uncertainty 
due to the power suppressed operator R 2 and the remaining //-dependence. 

One clearly sees that the change of the basis resulted in a considerable reduction 
of the theoretical error , almost a factor 3 in the case of Ar s /AM S ! 
To improve our theoretical knowledge of the mixing quantities further one needs more 
precise values of the non-perturbative parameters, like the decay constants or the 
power suppressed operators. If accurate non-perturbative parameters are available 
one might think about NNLO calculations (a 2 s or a^/mb-corrections) to reduce the 
remaining /i-dependence. 
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Figure 1: Error budget for the determination of Ar s in the old and the new basis. 



4.3 New Physics 

New physics (see e.g. [36]) is expected to have almost no impact on Ti 2 but it can 
change Mi 2 considerably. Therefore one can write 

r 12iS = rt 2 M s , M 12 , s = i\4 M s -A s ; A s = |A s |e^ (30) 

With this parameterisation the physical mixing parameters can be written as 



AM, 


= 2|M? 2 M S | 


•|A S | 






(31) 


Ar s 


= 2|ri2 )S | 


cos u 




*) 


(32) 


Ar s 


| r* 12 s 1 


COS 




^) 


(33) 


AM, 


\M®i\ 




[A.] 






Tl2 s | 


sin u 


b* M + < 




(34) 


|Mf 2 M s | 




|A S | 





Now we combine the current experimental knowledge about the mixing parameters 
to find out whether 5 s -mixing is described by the standard model alone, or whether 
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we already get some signals of new physics contributions 
The mass difference AM S is now known very precisely 



AM, = 17.77 ±0.10 {syst) ± 0.07 (^ps" 1 CDF. (35) 

For the remaining mixing parameters in the _B s -system only experimental bounds 
are available. The width difference Ar s /r s was investigated at ALEPH, BELLE, 
CDF, DO by analyzing the decays B s -> + [37], B s -> J/tf + [38] and 

Z? s — > K + + K~ [39] and the flavor specific lifetime of the -B s meson [10] • A recent 
combination of all these results yields [H] 



Ar s = 0.097 ± 0.042 ps" 1 , (36) 

Ar 

= (56 ± 24) x 10~ 4 . (37) 



AM, 



Except the angular analysis B s — > J/^> + <fi all other determinations are affected by 
some drawbacks, described in |7J. Moreover the DO collaboration has updated their 
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results [38] for the decay B s -> J/V + in (321 EH I2E! using lfb _ l of data. Setting 
the value of the mixing phase <p s to zero they obtain [121 [T71 [2S] 

Ar s = 0.12 ± 0.08^ ps" 1 , (38) 

allowing for a non-zero value of the mixing phase <p s they get 

Ar s = 0.17 ± 0.09 ± 0.03 ps" 1 , (39) 
S = -0.79 ±0.56 ±0.01. (40) 

In the following we will for Ar s and cj) s only use the numbers from Eq. ( |39i) and Eq. 
(HQ]). 

The semileptonic CP asymmetry in the B s system has been determined directly 
in [43] and found to be 

a s f irect = (24.5 ± 19.3 ± 3.5) • 10~ 3 . (41) 

Moreover the semileptonic CP asymmetry can be extracted from the same sign 
dimuon asymmetry that was measured in [H] to be 



a sl = (-2.8 ± 1.3 ±0.9) • 10" 3 . (42) 

Updating the numbers in (321 HE] one sees that 

a sl = (0.582 ± 0.030) a% + (0.418 ± 0.047) a s sl (43) 

In [l5l H6] the experimental bound for a d sl was used to extract from Eq. (142!) and 
Eq. (|43!) a bound on a s sl . Due to the huge experimental uncertainties in a d sl this strategy 
resulted in a large error on a s sl . Since in the ^-system there is not much room left 
for new physics contributions, we think it is justified to use the theoretical number 
of a d sl . Using a d sl = - (0.48 ± 0.12) ■ 10~ 3 we get from Eq.(g2J, Eq.gSJ and Eq.(|42D 
already a nice bound 



,dimuon 



;-6.0±3.2±2.2) • 10~ 3 . (44) 



Combining this number with the direct determination [13] we get our final experi- 
mental number for the semileptonic CP asymmetries 

a s sl = (-5.2 ±3.2 ±2.2) ■ 10~ 3 . (45) 

Now we combine these experimental numbers with the theoretical errors to extract 
bounds in the imaginary A s -plane by the use of Eqs. ( J3TT) . (|32l) . ( J33l) and (|34l) . see 
Fig. 0. 

The comparison of experiment and standard model expectation for AM S , AT S , 
<f) s , AT S /AM S and a s sl presented in figure [3] already shows some hints for deviations 
from the standard model. 
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Figure 3: Current experimental bounds in the complex A s -plane. The bound from 
AM S is given by the red (dark-grey) ring around the origin. The bound from 
AT S /AM S is given by the yellow (light-grey) region and the bound from aj s is given 
by the light-blue (grey) region. The angle (f)f can be extracted from AT S (solid lines) 
with a four fold ambiguity - one bound coincides with the x-axis! - or from the an- 
gular analysis in B s — ► J /^<p (dashed line). If the standard model is valid all bounds 
should coincide in the point (1,0). The current experimental situation shows a small 
deviation, which might become significant, if the experimental uncertainties in Ar s , 
a s sl and S will go down in near future. 
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5 Conclusion and outlook 



Theoretical predictions of the lifetimes of heavy mesons are in excellent agreement 
with the experimental numbers. We do not see any signal of possible duality viola- 
tions. To become more quantitative in the prediction of r#+ / T# d the non-perturbative 
estimates of the bag parameters B\, B 2 and ei, €2 have to be improved. For TB s /TB d 
more precise experimental numbers are needed, while in the case of tb c theoretical 
progress is mandatory. 

The theoretical uncertainty in the mixing parameter AM is completely dominated 
by the decay constant. We have presented a method (see [7] for more details) to re- 
duce the theoretical error in AT, AT/ AM and a,f s considerably. This relatively clean 
standard model predictions can be used to look for new physics effects in S s -mixing. 
From the currently available experimental bounds on Ar s and a/ s one already gets 
some hints for deviations from the standard model. This situation will improve dra- 
matically as soon as more data are available. 

For a further reduction of the theoretical uncertainty in Ar a much higher ac- 
curacy than currently available on the decay constants is necessary. If this problem 
is solved or if one looks at quantities like AT/ AM and aj s where the dependence 
on fs cancels, then the dominant uncertainty comes from the unknown matrix ele- 
ments of the power suppressed operators. Here any non-perturbative estimate would 
be very desirable. If accurate non-perturbative parameters are available one might 
think about NNLO calculations (a 2 s or a s /m;,-corrections) to reduce the remaining 
/^-dependence. 

I would like to thank the organizers of HQL2006 for the invitation and Uli Nierste 
for the pleasant collaboration. 
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